Tailoring Josephson coupling through superconductivity-induced nonequilibrium 
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The distinctive quasiparticle distribution existing under nonequilibrium in a superconductor- 
insulator-normal metal-insulator-superconductor (SINIS) mesoscopic line is proposed as a novel 
tool to control the supercurrent intensity in a long Josephson weak link. We present a description 
of this system in the framework of the diffusive-limit quasiclassical Green-function theory and take 
into account the effects of inelastic scattering with arbitrary strength. Supercurrent enhancement 
and suppression, including a marked transition to a 7r-junction are striking features leading to a 
fully tunable structure. 
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Nonequilibrium effects in mesoscopic superconduct- 
ing circuits have been receiving a rekindled attention 
during the last few years 0. The art of control- 
ling Josephson coupling in superconductor-normal metal- 
superconductor (SNS) weak links is at present in the 
spotlight: a recent breakthrough in mesoscopic super- 
conductivity is indeed represented by the SNS transis- 
tor, where supercurrent suppression as well as its sign 
reversal (7r-transition) were demonstrated 0, 0. This 
was achieved by driving the quasiparticle distribution in 
the weak link far from equilibrium 0. 0- through ex- 
ternal voltage terminals, viz. normal reservoirs. Such a 
behavior relies on the two-step shape of the quasiparti- 
cle nonequilibrium distribution, typical of diffusive meso- 
scopic wires and experimentally observed by Pothier and 
coworkers @. 

The purpose of this paper is to demonstrate that it is 
possible to tailor the quasiparticle distribution through 
superconductivity-induced nonequilibrium in order to im- 
plement a unique class of superconducting transistors. 
This can be achieved when mesoscopic control lines are 
connected to superconducting reservoirs through tunnel 
barriers (I), realizing a SINIS channel. The peculiar 
quasiparticle distribution in the N region, originating 
from biasing the S terminals, allows one to access several 
regimes, from supercurrent enhancement with respect to 
equilibrium to a large amplitude of the 7r-transition pass- 
ing through a steep supercurrent suppression. These fea- 
tures are accompanied by a large current gain (up to some 
10 5 in the region of larger input impedance) and reduced 
dissipation. The ultimate operating frequencies available 
open the way to the exploitation of this scheme for the 
implementation of ultrafast current amplifiers. 

The investigated mesoscopic structure (see Fig.QJ con- 
sists of a long diffusive weak link of length L j much larger 
than the superconducting coherence length (£q) oriented 
along the x direction. This defines the SNS junction of 
cross-section Aj. The superconducting terminals belong- 
ing to the SNS junction, labelled Sj (3 and 4), are kept at 




FIG. 1: Scheme of the Josephson transistor. The supercurrent 
I j (along the white dashed line) is tuned by applying a bias 
Vc across the SINIS symmetric line connected to the center 
of the weak link. All normal wires are assumed quasi-one- 
dimensional. 



zero potential. The SINIS control line is oriented along 
the y direction and consists of a normal wire, of length Lc 
and cross-section Ac , connected through identical tunnel 
junctions of resistance TZt to two superconducting reser- 
voirs Sc (1 and 2), biased at opposite voltages ±Vc/2. 
The superconducting gaps of Sj and Sc (Aj and Ac) 
are in general different. 

The supercurrent Ij flowing across the SNS junction 
is given by 0, @ 
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(1) 

and depends on the quasiparticle distribution function 
f(E). In Eq. Q, a is the normal-state conductivity 
which determines the normal-state resistance of the junc- 
tion according to Rn = Lj/crAj. The distribution 
function / reduces to the equilibrium Fermi distribution 
when Vn = 0. The energy-dependent spectral super- 
current @, 0, Imf^gj, can be calculated by solving the 
Usadel equations jlCj . Following the parametrization of 
the Green functions given in Ref. 0, these equations in 
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the N region can be written 

j E = -smh 2 (d)d x x, 
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d x jE = 0, 



HD ( 2 ) 
HDdlO + 2iEsinh0 + —(d xX ) 2 sinh(2(9) = 0, 

where D is the diffusion coefficient and E is the en- 
ergy relative to the chemical potential in Sj. 9(x,E) 
and \{x,E) are in general complex functions. For per- 
fectly transmissive contacts, the boundary conditions at 
the SjN interfaces reduce to 9 = arctanh( A j/E) and 
X = ±</>/2 in the reservoirs Sj, where (j) is the phase 
difference between the superconductors. 

As required by Eq. we must determine the actual 
quasiparticle distribution in the N region of the SINIS 
structure. This is controlled by voltage (Vc) and tem- 
perature, and by the amount of inelastic scattering in 
the control line. In the case of a short control wire with 
no inelastic interactions, the quasiparticle distribution, 
according to Ref. [TJ], is given by 



f(E,V c ) = 



Afl+J\f 2 '' 



(3) 



where A/i )2 = Ns c ( E ± eVfc/2) and Ti, 2 = T°{E ± 
eVc/2). The former are the BCS densities of states in 
the reservoirs Sc (labeled by 1 and 2 in Fig. QJ. T a {E) 
is the Fermi function at lattice temperature T [l2|. In 
this case Eqs. Q and yield the dimensionless transis- 
tor output characteristics shown in Fig. Ufa). The latter 
plots the supercurrent Ij versus control bias Vc at dif- 
ferent temperatures for a long junction (i.e., Aj 3> Ej, h , 
where E^ h — hD/L 2 j is the Thouless energy of the SNS 
junction, as this is the limit where the supercurrent spec- 
trum varies strongly with energy). We assumed <fi = n/2, 

t c c /tJ 



0.2, where T C C(J) are the critical tempera- 



tures of the superconductors $c(J) an d Lj such that 
Aj/E£ h = 300. 

At the lowest temperatures, by increasing Vc curves 
display a large supercurrent enhancement with respect 
to equilibrium at bias V c = 2 A c (T)/e = V C (T) (region 
I in the figure). Further increase of bias leads to a ir- 
transition (region II) and finally to a decay for larger 
voltages. This behavior is explained in Figs. 2(b,c,d) 
where the spectral supercurrent (solid line) is plotted to- 
gether with f(—E) — f(E) (dash-dotted line) for values 
of Vc and T corresponding to regions I, II and III, re- 
spectively. Hatched areas represent the integral of their 
product, i.e., the supercurrent Ij of Eq. (JIJ. In partic- 
ular, region I corresponds to the cooling regime where 
hot quasiparticles are extracted from the normal metal 
[Til ll^j . The origin of the 7r-transition in region II is 
illustrated by Fig. |5fc), where the negative contribution 
to the integral is shown. We remark that the intensity of 
the supercurrent inversion is very significant. It reaches 
about 60% of the maximum value of Ij at Vc — Vq (T) in 
the whole temperature range, nearly doubling the 7r-state 
value of the supercurrent as compared to the case of an 
all-normal control channel . In the high-temperature 
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FIG. 2: (a) Supercurrent vs control voltage Vc at different 
temperatures (T) for = tt/2 and T° = 0.2T C J (see text). 
Bias regions labeled (I), (II) and (III) indicate supercurrent 
enhancement due to quasiparticle cooling, high-voltage tt- 
state and low-voltage 7r-state in the high-temperature regime, 
respectively. These are qualitatively explained in (b), (c) and 
(d) where hatched areas represent the contribution to super- 
current arising in such bias ranges (see text). 



regime (T/T^ > 0.6), when the equilibrium critical cur- 
rent is vanishing, the supercurrent first undergoes a low- 
bias 7r-transition (region III in the figure), then enters 
regions I and II. This recover of the supercurrent from 
vanishingly small values at equilibrium is again the con- 
sequence of the peculiar shape of / (see Fig. EJcI)). No- 
tably, the supercurrent enhancement around Vq(T) re- 
mains pronounced even at the highest temperatures, so 
that I j attains values largely exceeding 50% of the junc- 
tion maximum supercurrent. This demonstrates the full 
tunability of the supercurrent through nonequilibrium ef- 
fects induced by the superconducting control lines. We 
remark that this is a unique feature stemming from the 
superconductivity-induced nonequilibrium population in 
the weak link. 

The length Lc of the SINIS control line can be addi- 
tionally varied to control the supercurrent by changing 
the effective strength of inelastic scattering in the N re- 
gion. For TZt ^> Rc, the distribution function f(E) in 
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the N region is essentially y-independent and we have 
1 



f(E))+N 2 (T2- f{E))] 



k J dujdeuj a X(uj,e, y,E) = 0. 



(4) 



Here vp is the normal-metal density of states at the Fermi 
energy, f2c is the volume of the N region and 1 is the 
net collision rate at energy E. At low temperatures, the 
most relevant scattering mechanism is electron-electron 
scattering [l^ and we can neglect the effect of electron- 
phonon scattering. Then 0, Il6| , 

e, y, E) = r>, s, y, E) - l° ut (u>, e, y, E), (5) 




2 eVr/A r 3 



and 



l m (cu, e, E) = (l- /( e ))(l - f(E))f(e - u)f(E + w),(6) 
e, £) = (1 - f(e - w))(l - /(£ + c))/(e)/(i?).(7) 

Electron-electron interaction is either due to direct 
Coulomb scattering 0, or mediated by magnetic 
impurities 15]. Below, we concentrate on the former, 
but the latter would yield a similar qualitative behav- 
ior. From the calculation of the screened Coulomb in- 
teraction in the diffusive channel, it follows 0] that 
a = —3/2 for a quasi-one dimensional wire and k — 
{■Ky/D/2h 3 / 2 v F A c )~ 1 0. We note that A c is the most 
relevant energy scale to describe the distribution func- 
tion for different voltages Vc- It is thus useful to replace 
lj — > lo/Ac and e — > e/ Ac, in order to obtain a dimen- 
sionless equation. Multiplying Eq. 10} by ^H-t^Icvf we 
obtain 



Mx{T x - }{E)) - Af 2 (f(E) - T$ = 
= K. coU I dudecj- 3 / 2 !^^^^^), 



(8) 



where K cM = f^^^/Ac" = V2f^j£-, R K = 

h/2e 2 and E! Th — HD/L C . In the absence of electron- 
electron interaction (JC C oU — 0) Eq. J3J) is recovered. 

The influence of inelastic scattering on Ij is shown in 
Fig. which displays the critical current of a long junc- 
tion at T — 0.1 Tf? for several values of /C co ;/. Here Ij 
is obtained by numerically solving Eq. (JSJ . The effect of 
electron-electron interaction is to strongly suppress the n- 
state and to widen the peak around V c . The 7r-transition 
vanishes for IC co u ~ 100, but the Ij enhancement due 
to quasiparticle cooling still persists in the limit of even 
larger inelastic scattering [13 • The disappearance of the 
7r-state can be understood by looking at the right in- 
set of Fig. which clearly shows how / (calculated at 
eVc — 2.5Ac) gradually relaxes from nonequilibrium to- 
wards a Fermi function upon increasing JC co a. The left 
inset shows how / (evaluated at eVc = 1.5 Ac) sharpens, 
thus enhancing Ij , by increasing JC co u . This effect follows 
from the fact that inelastic interactions redistribute the 
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FIG. 3: Supercurrent vs Vc for various KL co u with T = 0.1 
and — 0.2 T C J . Insets show the distribution function at 
eVc = 1.5A C (left) and eVc = 2.5A C (right) calculated for 
the same 1C co u values. 



occupation of quasiparticle levels in the N region, thus 
increasing the occupation at higher energy. As a con- 
sequence, higher-energy excitations are more effectively 
removed by tunneling, even for biases well below and not 
only around V c (as in the case of K. co u = 0). At the 
same time, supercurrent recovery at high temperature 
is gradually weakened upon enhancing K, co u. Notably, 
these calculations show that a rather large amount of in- 
elastic scattering is necessary to weaken and completely 
suppress the 7r-state. For example, using AI/AI2O3/CU 
as materials composing the SINIS line, JC co a — 1 corre- 
sponds to use a fairly long control line with Lc — 2.3 fj,m 

Changing the ratio T ( r shifts the Ij response along 
the Vc axis, the shape of the characteristics being virtu- 
ally independent of T c c . This translates into a differ- 
ent magnitude of control voltages Vc and power dissipa- 
tion P — IcVc, where lc is the control current across 
the SINIS channel. The function P(V C ) is plotted in 
Fig. Efa) for some ratios Tf/T/ at T = 0.01 T C J , assum- 
ing Tl T = 10 3 tt and T/ = 9.26 K. The impact of A c in 
controlling power dissipation is easily recognized. These 
effects clearly indicate that Ac <C Aj is the condition 
to be fulfilled in order to minimize P. In practice, the 
power dissipation for Vc > V c constitutes an experimen- 
tal problem as this energy needs to be carried out from 
the reservoirs. In a similar way the noise properties of 
the system are sensitive to the different TjF/T/ ratios. 
Assuming that the noise through one junction is essen- 
tially uncorrelated from the noise through the other, it 
follows that the input noise power S in the control line 
can be expressed as 

S = W~ I d£M[/(£)(l- ^i)+^i(l -/(£))]• (9) 
S(Vc) from Eq. 10 is shown in Fig. Efb) for the same 
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FIG. 4: (a) Power dissipated in the SINIS line vs Vc calcu- 
lated for various ratios T C C /T C J and T = 0.01T/. (b) Noise 
power S vs Vc calculated for the same parameters as in (a), 
(c) Differential current gain Gi vs Vc for T C C /T C J = 0.2. The 
inset shows Gj in the high-bias region. In all these calcula- 
tions we set K. co u = and T C J = 9.26 K (Nb). 

parameters of Fig. ^a). For example, for /T C J = 0.1 
(corresponding roughly to the combination Al/Nb), P 
obtains values of the order of a few 10 -15 W and S of 
some 10 -30 A 2 Hz _1 in the cooling regime, while these 
values are enhanced respectively to few tens of 10~ 12 W 
and 10 -26 A 2 Hz _1 for biases around the 7r-transition. 



In light of the possible use of this operational prin- 
ciple for device implementation, let us comment on the 
available gain and switching times. Input and output 
voltages are of the order of Ap/e and E^ h /e, respec- 
tively, so that it seems hard to achieve voltage gain. On 
the other hand, differential current gain G/ = dlj/dlc 
can be very large. For Vc > Vq a simple estimate gives 
Gj ~ (E^/A^CR-t/Rn), meaning that with realistic 
ratios H t /Rn (~ 10 3 ), G/ can exceed 10 2 . Gj(Vb) cal- 
culated for Tf/T/ = 0.2 is plotted in Fig. He) (the inset 
shows the gain in the 7r-state region). This calculation 
reveals that Gr can reach huge values, with some 10 5 
for Vc < Vq [21| and several 10 2 in the opposite regime. 
Remarkably, gain is almost unchanged also in the pres- 
ence of weak inelastic scattering (i.e., JC C oii = 1). The 
same holds for P and S. The highest operating frequency 
v of the transistor is limited by the smallest energy in 
the system: v < min ^{A c , Aj, E^ h ,E^ h , ^TZtC)^ 1 } 
where C is the tunnel junction capacitance. For an op- 
timized device, working frequencies of the order of 10 11 
Hz can be experimentally achieved in the high-voltage 
regime Vc >V C . For Vc < Vq, conversely, the response 
is slower (somewhat below 10 9 Hz), owing to the long 
discharging time through the junctions. 

We thank M. H. Devoret, K. K. Likharev, F. Pierre, 
L. Roschier, A. M. Savin, and V. Semenov for helpful 
discussions. This work was supported in part by MIUR 
under the FIRB project RBNE01FSWY and by the EU 
(RTN-Nanoscale Dynamics) . 



[1] See, for example, Theory of Nonequilibrium Superconduc- 
tivity, N.B. Kopnin (Clarendon Press, Oxford, 2001). 
[2] J. J. A. Baselmans, A.F. Morpurgo, B.J. van Wees, and 

T.M. Klapwijk, Nature (London) 397, 43 (1999); J. 

Huang, F. Pierre, T.T Heikkila, F.K. Wilhelm, and N.O. 

Birge, Phys. Rev. B 66, 020507 (2002); R. Shaikhaidarov, 

A.F. Volkov, H. Takayanagi, V.T. Petrashov, and P. Dels- 

ing, ibid. 62, R14649 (2000). 
[3] J. J. A. Baselmans, T.T. Heikkila, B.J. van Wees, and 

T.M. Klapwijk, Phys. Rev. Lett. 89, 207002 (2002). 
[4] A. F. Volkov, Phys. Rev. Lett. 74, 4730 (1995). 
[5] F.K. Wilhelm, G. Schon, and A.D. Zaikin, Phys. Rev. 

Lett. 81, 1682 (1998). 
[6] S.-K. Yip, Phys. Rev. B 58, 5803 (1998). 
[7] H. Pothier, S. Gueron, N.O. Birge, D. Esteve, and M.H. 

Devoret, Phys. Rev. Lett. 79, 3490 (1997). 
[8] See W. Belzig, F.K. Wilhelm, C. Bruder, G. Schon, and 

A.D. Zaikin, Superlattices Microstruct. 25, 1251 (1999) 

and references therein. 
[9] T.T. Heikkila, J. Sarkka, and F.K. Wilhelm, Phys. Rev. 

B 66, 184513 (2002). 
[10] K.D. Usadel, Phys. Rev. Lett 25, 507 (1970). 
[11] DR. Heslinga, and T.M. Klapwijk, Phys. Rev. B 47, 

5157 (1993). 

[12] At low control voltages there is a region of energies where 



A/"i = A/2 = 0. We assume that there the (otherwise weak) 
coupling to phonons makes the distributions at those 
energies equal to the equilibrium Fermi distribution. A 
more detailed discussion and another type of coupling is 
given in |T^|. 

[13] JR. Pekola, T.T. Heikkila, A.M. Savin, J.T. Flyktman, 

F. Giazotto, and F.W.J. Hekking, submitted. 
[14] M.M. Leivo, JR. Pekola, and D.V. Averin, Appl. Phys. 

Lett. 68, 1996 (1996). 
[15] A. Anthore, F. Pierre, H. Pothier, and D. Esteve, Phys. 

Rev. Lett. 90, 076806 (2003). 
[16] K.E. Nagaev, Phys. Rev. B 52, 4740 (1995). 
[17] B.L. Altshuler and A.G. Aronov, Zh. Eksp. Teor. Fiz. 75, 

1610 (1978) [Sov. Phys. JETP 48, 812 (1982)]. 
[18] A. Kamenev and A. Andreev, Phys. Rev. B 60, 2218 

(1999). 

[19] F. Giazotto, F. Taddei, T.T. Heikkila, R. Fazio, and F. 

Beltram, Appl. Phys. Lett. 83, 2877 (2003). 
[20] This is straightforward assuming as typical parameters 

TZt = 10 3 ft, D = 0.02 m 2 /s and A c = 200 ,ueV. 
[21] In this calculation we chose to include depairing by a phe - 

nomenological, but realistic parameter T = 10 _4 Ac uM- 

Its omission would lead to extremely higher Gi values. 



